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ABSTRACT 



Aims. We examine some properties of stars evolving close to the classical Eddington limit for electron-scattering opacity, 
when these stars maintain a chemically homogeneous structure as a result of mixing and/or mass loss. 
Methods. We consider analytical relations and models computed with the Geneva code. 

Results. Homologous, chemically homogeneous stars evolving with a constant Eddington factor obey a relation of the 
form p, 2 M = const. This applies, for example, to Wolf-Rayet (WR) stars in stages without hydrogen. The value of 
the constant may depend on the metallicity, initial mass, evolutionary stage, and physical processes included in the 
considered homologous evolutionary sequence. An average value of the constant between 20 and 40 in solar units is 
consistent with the masses of Galactic WR stars. 
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1. Introduction 

Stars must have a luminosity inferior to the Eddington 
limit, otherwise "... the radiation obser ved to be em i tted . .. 
would blow up the star" as written by Eddington (1926). 
This limit, which occurs when the outward acceleration of 
radiation pressure becomes equal to gravity, usually has a 
simple analytical expression (see Eq. [7]). Here, we present 
another simple expression of the Eddington limit applica- 
ble to a sequence of homologous stellar models that evolve 
without losing their homogeneous chemical structure. 

During their evolution, stars generally become chemi- 
cally inhomogeneous, hence our proposed expression would 
not be applicable. However, in some relatively rare cases, 
stars can indeed maintain a chemically homogeneous struc- 
ture during the H- and He-burning phases. This may result 
in massive stars underg oing a heavy m ass loss that removes 
the oute r stellar layers (Maeder 1980), or in strong internal 
mixing (|Maederlll987ri . 

The latest grids of stellar mode ls developed by the 
Geneva group ( Ekstrom et al.l I2012T ) show homogeneous 



evolutionary tracks for stellar masses above 50 Mq with 
even moderate initial rotation at 40% of the critical ve- 
locity. This confirms that depending on rotational mixing a 
bifurcation may occur in the general scheme of stellar evolu- 
tion: tracks evolve toward the blue instead of the usual way 
toward the red. For more rapid rotation, homogeneous evo- 
lution appears above lower mass limits. The homogeneous 
st ars rapidly enter th e Wolf -Ray et (WR) stages as s hown 
bv lMevnet fc Maeder! (|2005l) and lGeorgy et all (|2012D . 

The present work on the Eddington limit of homoge- 
neous stars is i n line and consistent with the results of 
iGrafener et al.l (|201l[ ) showing that " the proximity to the 
Eddington limit is the physical reason for the onset of WR 



type mass loss" . The expression we propose here may be 
particularly useful for deriving the mass loss rates of WR 
stars. 

An additional motivation for this study is that the ho- 
mogeneous evolution of massive hot stars has been pro- 
posed as the scenario l eading to long-soft gamma ray bursts 
(jYoon fc Langerl l2005). In this case, homogeneity would be 
due to the rapid axial rotation of the star. Homogeneous 
stars can also be formed during stellar mergin g. This sce- 
nario was reexamined bv lGlebbeek et all (|2009f ) as a possi- 
ble means of forming intermediate-mass black holes in glob- 
ular clusters. 

In Sect. [5J we proceed to some basic analytical develop- 
ments. In Sect. [3j we consider a polytropic structure and 
in Sect. H]the case of homologous evolution. WR observa- 
tions are examined in Sect. [5l while Sect. [5] presents our 
conclusions. 



2. The upper mass limit for homogeneous stars 

2.1. The mass-luminosity relation 

We begin with and develop briefly some basic relations. The 
equation of radiative equilibrium 



dT 4 



■iirr 2 Sup dr ' 
implies the scaling for the luminosity L 

up 



(1) 



(2) 
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where k and p are the mean values of the opacity and den- 
sity in a star of radius R and central temperature T c . The 
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Fig. 1. Mass-luminosity re lation for ZAMS stars with X = 
0.70 and Z = 0.014 (from iGrafener et al.|[20TTl in black), 
compared to the limits given by Eq. [5] for the same compo- 
sition and various T-values as indicated. 



equation of hydrostatic equilibrium 4r- = —pg gives a scal- 
ing for the central pressure 



Pc (X 



R* 



(3) 



The equation of state for massive hot stars is that of 
a mixture of perfect gas and radiation P = P gas /f3 = 
kpT/(f3pm u ), where j3 is the ratio of the gas pressure to 
the total pressure. This implies that 



rp M/?C p 

l c (X P c oc 



(Pc/p) R ' 



(4) 



As the star is chemically homogeneous, the mean molecular 
weight p is assumed to be constant throughout, (for cool 
stars p may vary owing to changes in the degree of ioniza- 
tion). For this expression of T c , the luminosity L behaves 
like 

~ 4 M 3 



Pc 



(5) 



where A is a proportionality constant. The above scalings 
are rather rough, thus A is almost constant only as long 
as the model structures do not differ too considerably. We 
consider below the polytropic and homologous structures. 
The term /3 and the ratio of the two densities are often 
omitted from the mass-luminosity relation. 



2.2. The Eddington limit 

According to the definition, the Eddington luminosity is 
given by 

1 dP rad re s Ledd GM 

(6) 



Ps 



dr 



4vri? 2 



R 2 



The first equality expresses the gradient of the radiation 
pressure in term of the radiative flux, while the second says 
that this gradient is equal to the modulus of gravity. The 



subscript "s" stands for quantities at the surface. This gives 
the classical Eddington limit 



A-kcGM 



^EDD 



(7) 



Here, it is the opacity re s in layers close to the surface which 
intervenes. However, in hot massive stars, the main opacity 
source is electron scattering. To first order, the opacity is 
the same throughout a homogeneous star, thus here we no 
longer distinguish the two re values. 

The maximum luminosity effectively observed Ly may 
lie below the formal limit Ledd and be given by 



L r = TL 



EDD 



(8) 



which is the effective limit for some Y < 1. 

We consider the possibility that the real Y is lower than 
1.0, since there are several physical reasons for that. Firstly, 
the actual opacities in the outer stellar layers are higher 
than just the electron-scattering opacities assumed for the 
classical Eddington limit. In particular, strong resonance 
spectral lines are driving the heavy radiative mass loss of 
massive stars. Any pulsational instabilities, that is either 
radial or non-radial oscillations, are favoured by the de- 
crease in the adiabatic exponents Y\, Y2 and Y3. These 
exponents are all reduced by an increase in the ratio of the 
radiation pressure to the gas pressure, as it occurs close to 
the Eddington limit. The excitation mechanisms of either 
p modes or g modes, as observable in asteroseismology, are 
also enhanced by a high radiation pressure. This is also the 
case for the strange modes develo ping in very massive star s 
with a high radiation pressure ()Papaloizou etaLl fl997h . 
Thus, there are several effects and instabilities that may 
place the real maximum luminosity of stars of a given mass 
and composition below the academic value of the Eddington 
limit. These theoretical considerations are in close agree- 
ment with observations. For example, all the most luminous 
WR stars in the Arches cluste r have T-values between 0.2 
and 0.45 (Gr afener et al.ll201.il) . i.e. at some significant dis- 
tance below the formal limit. 

We emphasise that the effective Y may depend on the 
parameters of the considered stellar population, such as its 
metallicity Z, etc. For example, for a higher metallicity Z, 
the actual opacity in the innermost layers of the wind is 
larger, implying that the momentum of radiation onto the 
outer stellar layers is larger, and also that pulsational in- 
stabilities due to the re-mechanism are enhanced. Thus, the 
effective Y is likely to be lower at higher Z. Other effects 
influencing the mass loss, such as rotation or tidal effects 
in binary system, may also influence the effective Y. 

2.3. The effective limit 

We assume that the highest stellar luminosity L given by 
the mass-luminosity (M — L) relation \E\ is equal to Ly de- 
fined by Eq. [SJ 



. 4a4 p A M 3 ^AttcGM 

Pc K 



Y 



(9) 



Figures [T] and [5] show the lines corresponding to various 
r and the mass-luminosity relation from numerical mod- 
els, respectively, for stars on the zero-age main sequence 
(ZAMS) and for WR stars without hydrogen. In both cases, 
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Fig. 2. Mass-luminosity relation for WR stars without hy- 
drogen (in black) compared to the limits given by Eq. [5] 
for the same composition and various T-values as in- 
dicated. The continu ous black line refers to models by 
iGrafener et al.l (|2011| ) with X = 0, while the broken black 
line represents the ave rage mass-luminosity of W R stars 
without hydrogen from ISchaerer fc Maederl (|l992f) . 



we see that for T = 1 the M — L relation becomes parallel to 
the limit for high masses, implying that the limit is never 
reached, or in other words is only reached for extremely 
high masses or a value of /? tending toward zero. Therefore, 
Eq.Ej]for r = 1 would not correspond to realistic cases. The 
factor r measures the proximity to the Eddington limit. For 
T-values lower than 1.0, the M — L relation may intersect 
a Lr-line for some realistic stellar masses Mr. For exam- 
ple, the line T = 0.25 intersects the M — L relation for a 
star of about 57 Mq on the ZAMS (Fig. [J), while for stars 
without hydrogen the same T-value is reached by a star of 
about 12 M (Fig.©. 

Relation ((9]) defines the maximum mass Mr allowed for 
a given value of T 



f Mr = 



1/2 



with A in cm 4 

, ,'A-kcG 
with 



1 

K 

s' 3 

1/2 



(10) 



m g 



for chemically homogeneous stars with electron-scattering 
opacity. The first parenthesis on the right-hand side is a 
constant. We see that in general the product /j, 2 M is limited 
by a product of terms, where the last two are not constant. 
This product is examined in the next few sections. 



3. The polytropic case n = 3 

We consider equation (fit)]) for some particular models. We 
first examine the case of a polytropic structure with a pres- 
sure of the form P = K g l+l / n and index n — 3. If internal 
variations in /3 could be neglected, this case would be a 
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Fig. 3. Variation of j3 as a function of stellar mass for mod- 
els on the ZAMS with X = 0.705 and Z = 0.02. The max- 
imum in outer layers and th e minimum va lue in the centre 
for a given mass are shown (jMacdcr 200^). The polytropic 
case given by Eq. [TJ]lies between the two curves. 



good approximation for stars with an equation of state of 
perfect gas and radiation, since the pressure P could be 
written 

i 

4 



P = 



k 



a \/j,m u 



1-P 

/3 4 



(11) 



For such a polytrope, the equilibrium equations shows that 
the ratio of a star of mass M to the maxim um mass M max 
for this type of objects is feddingtonl ll926D 



GM 



(4 if) 3 
AttG 



(12) 



With the expression of K given by the bracket term with 
power 1/3 in Eq. (fTTj) . this becomes 



1-/3 = 



3 • 4 3 fc 4 M2 . 



(13) 



The parenthesis contains only constant terms and is equal 
to 7.562 ■ IP" 70 g~ 2 where we have taken M max = 2.018 
(|Chandrasekharlll939h . 



Expressing the mass in solar units, we find that 
(1-/3)V2 



/3 2 



0.0547 \i (M/M ) 



(14) 



This variation in the parameter /3 as a function of M is 
illustrated in Fig. [31 where it is also compared to the central 
and external values of realistic stellar models on the ZAMS. 
The above relation indicates that, in a polytrope, the ratio 
j3 of the gas to the total pressure is a function only of n 2 M. 
This function is illustrated in Fig. [4] 

In a polytrope with a given index n, the ratio of the cen- 
tral to the average density is always the same. For n = 3, 
the ratio (p c /p) is equal to 54.2. Thus, besides the term 
r 1 / 2 , the right-hand side of Eq. [TU] contains only constant 



3 



Maeder et al.: The Eddington limit 



terms and a function /? = f3(fi 2 M). This function is mono- 
tonically decreasing (Fig. [4]). For high masses (say above 
80 Mq), the parameter ft tends to vary as (^ 2 M) -1 / 2 . 
Thus, the right-hand side of Eq. [TUJ which contains the 
term (1//3) 2 , behaves for very high masses like M . Thus, 
Eq. [TU] has no solution in this limit. This is the same lim- 
iting situation as for the usual Eddington limit with T = 1 
illustrated in Figs. Q] and [5J 

However, for not too large values of p 2 M, the function 
f3~ 2 , which appears on the right-hand side of Eq. 1101 grows 
faster than for very high masses. Thus, for some ranges of 
r -values smaller than 1.0, the curves representing the left 
and right-hand sides of Eq . [TOl may cross each other. This 
defines the maximum possible mass Mr for the considered 
F- values. 

The solution also depends on the value of the param- 
eter A in the mass-luminosity Eq. [5j This relation is not 
defined by the polytropic structure, which just results from 
the mechanical equilibrium expressed by Poisson's equa- 
tion. Nevertheless, the polytropic case is extremely interest- 
ing, since it shows that the product /i 2 M may be bounded. 




Fig. 4. The relation expressing f3 as a function of /i 2 M for 
a polytrope with index n = 3, as given by Eq. 1141 where M 
is expressed in solar units. 



4. Homology relations 

The polytropic case is rather restrictive and we now con- 
sider the more general case of an homologous sequence of 
chemically homogeneous models. These are models in which 
some parameters, such as the density and linear dimensions, 
are multiplied at each point by some con stant factors to 
produ ce another equilibrium configuration (jChandrasekharl 
1939). This typically correspond to a sequence of models 
with similar structures that differ by a slowly varying pa- 
rameter. Here, the changing parameter is the mean molec- 
ular weight fj, in an evolutionary sequence of chemically ho- 
mogeneous stars. We examine the implications of the limit 
given by Eq. (|TU| for such a sequence. 

First, in a homologous sequence of stellar models with 
the same density profile, the ratio (p c / p) remains constant. 



We now have to examine the behaviour of /3. We make a ho- 
mology transformation r = constant and M[ — xM r . The 
equation of hydrostatic equilibrium dP/dr = — {GM r /r 2 ) g 
implies that 

P' =x 2 P . (15) 

Thus, in Eq. Qj] on the left-hand side P behaves like x 2 
while on the right-hand side, the density term g 4 / 3 behaves 
like a; 4 / 3 . To be consi stent, the square bra cket in Eq. [TTI 
must behave like x 2 / 3 (jChandrasekharl I1939D . i.e. 



1-/3 

/? 4 



M 2/3_ 



This means that 



(1-/?) 1/2 
/3 2 



li 2 M 



(16) 



(17) 



Thus, we have returned to an equation similar to Eq. [T4l 
which applies to any point in the star, including the centre. 

The general dependence of j3 on /i 2 M can also be verified 
by expressing the ratio of the radiation to gas pressures 



Pr 



ad 



gas 



1 a T 3 pm u 
3 g k 



(18) 



We take the average density and average temperature T = 
§ ii ir L 1 TF'" wnere 9 is a concentration factor, with q = 
3/5 for an homogeneous density and q = 3/2 for a polytrope 
with n = 3. The above ratio becomes 



-^rad 

' (>';is 



f a(^)% 3 G WM 2 



(19) 



Since P ra d/-P g a S = (l-p)//3, one gets finally (|Maederi r2009') 



2 



(1-/3)V2 



Ana 



1/2 



k 



1 1 

G 3 / 2 q 3 / 2 



(20) 



There are only constant terms in the second member. For 
q = 3/5 (homogeneous density), the second member is 
equal to 21.93 Mq and for q — 3/2 (polytrope n = 3), 
it is equal to 5.55 Mq. Thus, we have verified in two ways 
that /3 is a function only of p?M. As mentioned above, this 
necessarily implies that for some range of T-values smaller 
than 1.0, the term p?M is bounded by a constant. 

We estimate this constant, for example, for a 120 Mq 
star on the ZAMS with X = .68, Z = 0.02, \x = 0.627, 
log(L/iQ) = 6.252, central and average densities of 1.483 
and 0.0438 g cm- 3 , and (3 = 0.495 (|Maederll200l . These 
parameters give on the left side of Eq. [TU]a value [i 2 M = 
47.18 Mq, while on the right side we have 



/4ttcG 



1/2 



— I I -r I = 75.98 Mr. 



© ■ 



(21) 



The equality of the two sides of Eq. [TTJ is realised for T = 
0.386. This value is consistent with the T-value of a 120 
Mq in Fig. \T\ This means that if this T- value were to be an 
effective limiting value, the limit fi 2 M — 47.18 Mq would 
not be overcome. For a pure helium star, the corresponding 
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Fig. 5. Mass- luminosity relations of evolutionary models 
with initial X = 0.72 and Z = 0.014 for massive stars with 
an initial rota tional velocity equal to 40% of the critical 
velocity ffrom lEkstrom et aDl2012h . 



limit Mr, for this particular T, would be 26.55 Mq instead 
of 120 Mq on the ZAMS. 

As another example, we see from Figs. [1] and [2] that, 
for r = 0.25 (where the intersection of the T-line with the 
M-L relation is well-defined), the ratio of the two masses 
57 Mq (with p = 0.617) and 12 Mq (with p = 1.340 ) 
corresponds to the inverse of their /z 2 -values. 

A crucial question now is whether during some phase the 
WR stars evolve keeping the same value of T. Fig. \b\ shows 
the mass-luminosity relation for rotating stars of differ- 
ent m asses evolving away from the ZAMS (jEkstrom et al.l 
l2012f) . The evolution is followed until the end of the cen- 
tral C-burning phase. The location of the Eddington limit 
for stars without hydrogen is indicated. We see that rather 
long parts of the evolutionary tracks in the phases of WR 
stars without hydrogen (types WNE and WC) are parallel 
to the Eddington limit at the corresponding masses. These 
tracks typically have a luminosity that is lower than the 
Eddington limit by a factor of about two. They correspond 
to a line T w 0.5, (the track for 85 Mq corresponds to 
T = 0.6, but is still parallel). 

Since these tracks have been calculated with mass-loss 
rates based on observations, the aforementioned parallelism 
indicates that the F-factors of WR stars remain constant 
during the concerned parts of the evolution. Incidentally, 
we note that the above value of T = 0.5 is an upper bound 
to the values T = 0.20 -0.45 observed in the Arches cluster 
(|Grafener et al.l [201 lh ■ This is consistent with the stellar 
models of Fig. [5] seeming to predict too high WR lumi- 
nosities with respect to the observed ones. This may be an 
indication that the mass-lo ss rates applied in the models 
are too low at some stages (jGeorgy et al.ll2012l) . 

As a conclusion of the last two sections, we note that, 
in a sequence of homologous chemically homogeneous WR 



models, the product p 2 M is effectively bounded by a con- 
stant factor 

M 

M 2 T7- < Cedd , (22) 

Mq 

where we call Cedd this constant. In the next section, we 
try to estimate the value of Cedd from observations of WR 
masses in the solar neighborhood. 

Here, it is very important to realise that Eq.lTUland con- 
sequently Eq. [52] are just logical consequences of the classi- 
cal Eddington limit for electron-scattering opacity. Thus, 
the above results do not contain subtle consequences of 
physical effects that are not present in the hypotheses. For 
example, if we considered such an effect as the metallic- 
ity dependence of the WR mass-loss rates and if this effect 
were to arise from bound-free and line opacities, it would 
be normal that our results do not contain it. Similarly, if 
some effects in the mass-loss rates originated from rotation, 
magnetic field, or binarity, it should also not be expected 
that the present results reproduce them. To summarise, our 
results are just equivalent forms of the classical Eddington 
limit based on electron-scattering opacity for some specific 
assumptions. In particular, Eqs.rT0]or[22]only apply to poly- 
tropes of a given n-index and to evolutionary sequences of 
chemically homogenous stars evolving homologously. For 
cases that significantly deviate from these hypotheses, our 
conclusions do not apply. 

Therefore, Cedd is not a "universal constant" for 
all stars. For example, evolutionary sequences of dif- 
ferent metallicity Z are likely to have different val- 
ues of Cedd for two reasons. Firstly, the product 
(4ttcG/A) 1/2 (pc/pf (1//3 c ) 2 may be slightly different. 
Secondly, the effective T expressing the ratio of the real ob- 
served maximum luminosity to the theoretical Eddington 
luminosity for electron-scattering opacity may be differ- 
ent as discussed in Sect. 2.3. The same kind of remark 
may apply to other differences of input parameters. Within 
homologous sequences of homogeneous models with dif- 
ferent input parameters, each sequence may evolve with 
p 2 M = constant, but with different values for the constant. 
However, polytropes with n = 3 all have the same constant 
as shown in Sect. 3, thus departures from this polytropic 
relation is an interesting point to check. 

Incidentally, we remark in Fig. [5] that the strong bright- 
ening of the WR stars in the very late stages (from the core 
contraction that precedes central C burning) brings the star 
very close to the Eddington limit. This suggests that the 
WR stars may experience extremely high mass-loss rates 
in the relatively fast stages preceding the supernova explo- 
sion. When this strong mass loss occurs, just before the 
SN explosion, it may hide from view the progenitors em- 
bedded in a d ense envelope. This may be the reason why 
ISmarttl ([2009T ) reports no detection of progenitors for ten 
Ibc supernovae having pre-explosive images. 



5. A quantitative estimate based on WR stars 

For homog eneous stars near the Eddington limit such as 
WR stars (|Grafener et al.l 1201 ll ). the condition in Eq. 
implies that the rate of mass loss dM/dt is determined by 
the growth rate of the mean molecular weight p resulting 
from the internal nuclear reactions 



dM 
~dT 



M dp 
p dt 



(23) 
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Interestingly enough, we see that the exact value of Cedd 
does not play a role in the amplitude of the mass-loss rates. 

However, the value of Cedd is still important to deter- 
mine at which stage of evolution the condition in Eg. Inland 
rates Eq. [23] start to apply. For a star on the ZAMS, the 
value of the product /i 2 M may initially be below the limit 
Cedd- However, as a result of nuclear reactions the value 
of [a may increase up to a point where the product n 2 M 
reaches the limiting Cedd- From there, in the sequence of 
homogeneous homologous models, the mass-loss rates will 
be controlled by Eq. [23] This does not necessarily mean 
that, for a given WR mass, the mass- loss rates will always 
be the same, since depending on Cedd the relation between 
M and fj, will not be the same and the rate of change d/i/dt 
may also be different. 

To estimate the effective constant Cedd, we now con- 
sider various groups of stars at different stages of evolution, 
studied in most recent literature. We first consider O-type 
stars on the ZAMS, which evidently are homogeneous with 
typically a ^t-value of about 0.62. On the basis of a large 
sample of 90 clusters, Ma schberger fc Clarkd (|2008f) deter- 
mined the masses of the most massive stars, to examine 
whether these are related to the cluster mass. The most 
massive star they found has a mass of 80 Mq. This is 
lower than the extreme cases of some LBV and H-rich WN 
stars, which have a lumi nosity corresponding to mass es of 
up to 120 Mq or above (jHumphrevs fc Stanekll200l . On 
the other hand, the lowest in itial stellar masses leading to 
WR stars are arou nd 25 Mq (|Maederlll98lHCrowtherll2007t 
ISander et aTll2011[) . Thus, for ZAMS stars we consider that 
the progenitors of WR stars lie in the range of from 25 Mq 
to about 120 Mq, most of them originating from stars with 
initial masses below 80 Mq. 

A most useful constraint comes from the WN stars with- 
out evidence of hydrogen at their surface. These stars are 
essentially pure helium stars with a mean molecular weight 
of 4/3, although some of them might already have both C 
and O deep in their interior that has not yet reached the 
surface. However, at this stage the products of He burning 
would relatively soon appear in the outer layers. Thus, it is 
likely that these WN stars have a structure close to chemi- 
cal homogeneity with an average [i-vahie not far from 4/3. 
A large sample of WN s tars without hydroge n has been 
analysed in the study bv lHamann et al.1 (|2006l) . We found 
that 30 of them have actual masses between 13 Mq and 
35 Mq with an average of 18.5 Mq. There is only one 
star with a mass above 25 Mq, thus most are in the range 
of from 13 Mq to 25 Mq. Masses of WN stars in binary 
systems were also determined. The lower bound is a value 
of 10 Mq, and the upper bound is not meaningful since 
large WN stars of large masses are generally hydrogen rich 
(|Crowtherll2007h . Thus, W N binaries have about the same 
lower mass as in the data of lHamann et all with a marginal 
trend to be slightly lower. 

Another interesting constraint is provided by WC stars. 
These are essentially stars in the process of helium burn- 
ing, that is to say they are helium stars with already a 
certain amount of C and a smaller amount of O. In view 
of the lifetimes of the nuclear burning phases, only about 
1% of WC stars may already be in the stage of central C 
burning and develop strong chemical inhomogeneity. Thus, 
most of the WC stars are likely to have a /x-value between 
that of ionised helium, i.e. 1.33 and that of ionised carbon, 
i.e. 1.71 (for ionised oxygen /i = 1.78). There may be some 
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Fig. 6. The curves yU 2 M = Cedd for 3 values of Cedd- 
The location of the ZAMS, He-stars and CO-cores are indi- 
cated. The average point for 30 W N stars without hydrogen 
studied bv lHamann et all (|2006fl is shown, with indications 
of the limits within they are observed. From the study by 
ISander et all (l20Tl . we show the averages for 4 transition 
WN/WC stars, for 39 WC stars (with indications of the 
limits) and for 2 WO stars. 



departures from a constant internal mean molecular weight 
in the WC phase. However, these are probably not large in 
view of the heavy mass loss and the relatively limited range 
of /i- values present. The stellar paramet ers of 39 WC stars 
were determined by ISander et al.l (|201lh . The average ac- 
tual mass, as determined from the luminosities, is 12.9 Mq. 
This is quite consistent with the values of the masses of 
WC stars determined fro m binary WR sta rs. These masses 
range from 9 to 16 Mq (|Crowtherl I2007D . There are also 
four transition WN/WC stars with an average mass of 21.5 
Mq and two WO stars with an avera ge mass of 22 Mq. 
As pointed out by ISander et al.l (|2011l) , the WO stars may 
stem from the most massive stars. Their location in Fig. [5] 
would agree with this conclusion. 

These various data are represented in Fig. [6j which 
shows masses versus mean molecular weights. We indicate 
the curves corresponding to Cedd= 20, 40, and 60. We see 
that the data based on significantly large samples are con- 
sistent with a general average value of Cedd between 20 
and 40, if we account for the weight of the samples of WN 
and WC stars. It is indeed uncertain whether the evolu- 
tionary sequence from O-stars to WN and WC stars forms 
a unique homologous sequence. This path may be "tan- 
gential" to one or two different homologous sequences with 
slightly different values of Cedd- The more the homogene- 
ity is satisfied, the more of course the homology property 
will be verified. This means that the observationally esti- 
mated value of Cedd between 20 to 40 is an average value 
covering the evolution from the MS of massive stars to WR 
stars and that in the future it might perhaps be more clearly 
specified for some parts of the evolution. 

Figure El also shows that the above observations are con- 
sistent with a curve of the type given by Eq. [22] although 
the WR stars originate from a rather broad interval of ini- 
tial masses. These results are in agreement with the mass- 



es 
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loss rates of WR stars being gove rned by the proximity to 
the Eddington limit, as stated bv lGrafener et all (]201l| ). 



6. Conclusions 

We have shown that, in a sequence of homologous, chem- 
ically homogeneous stars, a relation of the form of Eq. [22] 
is followed, if the star keeps a constant T-factor. Different 
homologous sequences characterised by different input pa- 
rameters may have different values of the constant in Eg. |2"2"1 
Recent determinations of WR masses of various subtypes 
tend to support an evolution according to the condition 
given in Eq. [22] and allow us to estimate that the aver- 
age value of the effective numerical constant in Eq. [22] is 
between 20 and 40. 

The validity of our aforementioned hypotheses (homo- 
geneity, homology, constant T) must of course be verified 
in the domain where the above relation would be applied 
in evolutionary calculations. 
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